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ABSTRACT 



A second-order perturbation theory approach is developed to calculate temperature 
anisotropies in the cosmic microwave background. Results are given for open universes 
and fluctuations corresponding to CDM models with either Harrison- Zeldovich (HZ) or 
Lyth-Stewart-Ratra-Peebles (LSRP) primordial energy-density fluctuation power spec- 
trum. Our perturbation theory approach provides a distinctive multipole contribution 
as compared to the primary one, the amplitude of the effect being very dependent on 
normalization. For low-0 models, the contribution of the secondary multipoles to the 
radiation power spectrum is negligible both for standard recombination and reionized sce- 
narios, with the 2-year COBE-DMR normalization. For a flat universe this contribution is 

0.1 — 10% depending on the reionization history of the universe and on the normalization 
of the power spectrum. 



Subject headings: cosmology: cosmic microwave background - O < 1 - 2nd order perturba- 
tions 
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/. INTRODUCTION 



In the absence of a cosmological term (A = 0) and assuming standard recombination, 
anisotropies in the temperature of the Cosmic Microwave Background (CMB) on large 
angular scales (> {20}^'^)°) are generated via the Sachs- Wolfe effect (f2 = 1), (Sachs and 
Wolfe 1967), or the generalized Sachs- Wolfe effect (fi < 1). This last case includes not 
only gravitational fluctuations at recombination but the integrated effect that depends 
on curvature of the time-varying potential along the photon trajectory (Anile and Motta 
1976, Abbott and Schaefer 1986, Traschen and Eardley 1986, Gouda et al.l991). 

For smaller angular scales, there are primary contributions coming from recombina- 
tion: photon fluctuations and the Doppler effect from last scattering. However, secondary 
anisotropies can also arise in the microwave sky due to the following physical processes: i) 
differential gravitational redshifts and blueshifts of the growing non-linear density fluctua- 
tions act on photons travelling to the observer that contribute an integrated gravitational 
effect (Martinez- Gonzalez, Sanz and Silk, 1990, 1992, 1994), ^ ii) hot gas inhomogeneously 
distributed between recombination and the present epoch generates anisotropy associated 
with non-linear flows (Vishniac, 1987), and iii) dust present in the early universe, produced 
in the process of formation of bound structures, also generates anisotropies especially in the 
Wien region (Bond et al., 1991). Such effects are important if the last scattering surface 
is modifled by reionization. An interesting possibility is that the secondary temperature 
fluctuations imprint signatures that differ in angular scale from the primary fluctuations 
(see section III). 

We will focus in this paper on the flrst of the above mentioned effects: the generation 
of anisotropies by 2nd order density perturbations in an open model. For the fluctuations, 
we will considerer a CDM model with null/negligible baryon content. Recent papers on 
this topic deal with the hypothesis of a flat background: 2nd order perturbations were 
considered by Martinez-Gonzalez, Sanz and Silk (1992) whereas some (analytic) models 
for the non- linear evolution were treated by Martinez-Gonzalez, Sanz and Silk (1994). 
Numerical simulations dealing with full non-linear evolution for (more realistic) models 
such as CDM have been performed by Tuluie and Laguna (1995). Seljak (1995) has 
also done numerical calculations based on N-body simulations performed by Gelb and 
Bertschinger (1994) for a flat universe. The generic conclusion one can extract from this 
work is that, at least for f2 = 1, the level of anisotropy in the CMB generated at low-z 
by quasi-linear and non-linear evolution of the matter is AT/T ^ 10~^, except for some 
unrealistic models where growth commences at very high z. 

We have computed the resulting fluctuations, which are of second order in perturbation 
theory relative to the uniform open cosmological background, for the case in which the 
density fluctuations are Gaussian. We expect non-linear contributions to further increase 

tft is helpful to distinguish between the integrated Sachs- Wolfe effect (Sachs and Wolfe 
1967; Hu and Sugiyama 1995) that occurs across linear fluctuations only when Q 1 from 
the differential effect across clusters (Rees and Sciama 1968) that is present even if Q = 1. 
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these minimal fluctuations, but such effects are highly model-dependent. Our perturbation 
theory approach provides a relatively robust prediction: for CDM open models, 5T/T has 
a distinctive multipolc contribution compared to the one expected in standard inflationary 
models at the linear level. The amplitude of the effect is very sensitive to normalization. 

//. THE SECOND ORDER GRAVITATIONAL EFFECT 

a) The integrated gravitational effect: 

We obtained (Martinez-Gonzalez, Sanz and Silk, 1990) an expression for the secondary 
anisotropics generated by the linear gravitational potential (fi{t,x) associated with non- 
linear density fluctuations A{t,x). For a flat or open universe with vanishing pressure, 

i^T/T)^^^^^^^^^ = 2j\t^{t,x) , vV = 61^a-^A(t,f) , (1) 

where the scale factor a{t) is normalized at the present time (ao = 1) and we choose units 
such that c ~ SnG — 1 and the Hubble length at the present time do = 2if~^ = 1. The 
previous equation giving the integrated gravitational effect is valid in the open case only for 
scales below the curvature scale. The line integral must be performed along the geodesic 
associated with the flat or open Friedmann background 

x{t,n)^X{a)n , a(A) ^ ^ J| ; (2) 

where n is the unit vector in the direction of observation. This integrated effect, except 
for the factor 2 that comes from general relativity, can be understood in terms of New- 
tonian mechanics: it represents the work performed by the photons, propagating from 
recombination to the present time, against the non-static gravitational potential (p. 

b) Calculations: 

Let us consider either an Einstein-de Sitter or open universe as background. Perturba- 
tion theory up to 2nd order, for vanishing pressure, gives the following expression for the 
density fluctuations: 

A{t,x) = D5^D^[^5^ + V5-Vi+'^-i,i^i^^^] , V^i = 5 (3) 

where 5{x) = 5^(1 + -^r)- The previous expression for A is exact up to second order for 
a flat model (Peebles, 1980), whereas it is a good approximation for open models with 
n>{).l (Bouchet et al. 1993; Catelan et al. 1995). 

Consequently, the second order effect, as given by equation (1), amounts to 
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(^^/^).econda.. = 12^^ da <^{x{a)){^)\2f - 1) 



where 

V20(f ) = ^ W • + ^O.-e^^' • (4) 



In the previous equation, D[a) is the growing mode of the perturbations normahzed to 1 
at the present time. For a flat universe: D = a, whereas for an open universe (Peebles, 
1980): 

n 9{x) A s dlnD 

9{x) = l + {l)[l + {l + l)'/Mi^ + xV/'-x'/')] . (5) 

Jb Jb 

The basic function to be calculated for any experimental set-up is the multipole com- 
ponent Ci =< af^ > (where temperature fluctuations expanded in spherical harmonics 
Yim are given by AT(n) = I]« Zlm=-« "«rn5^«m(^)), given by the following expression for 
the second order gravitational effect (we shall not consider the monopole and dipole in the 
calculations of the temperature anisotropy below): 



Here, p = A/[l — (1— f2)A^], ji is the Bessel function of fractional order and A^- is the distance 
from the observer to the recombination surface. The function P(2) (k) is the power spectrum 
associated with the 2nd order density perturbation 52 = + V5 ■ + and is 

related to the power spectrum of the time derivative of the potential ^ by = ■^P{2) (k)- 
The second order perturbation power spectrum P(2)(^) is given in terms of the first order 
power spectrum P{k) by the equation (Goroff et al. 1986, Suto and Sasaki, 1991) 



In the limit of small k, there is a cancellation of the three terms contributing to 62, implying 
that P{2)ik) has very little power on large scales. In particular, the generic behaviour is 
P(2){k) oc k'^ for small A;, independently of the primordial power spectrum. 

With regard to a possible secondary contribution coming from 3rd order density per- 
turbations for open models (in the case of flat models this does not exist because the 1st 
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order gravitational potential is static), we have the following comment: the coupling of the 
1st order gravitational potential ^ with the 3rd order potential ip^^^ oc ^ gives a 

kernel for the integrated gravitational effect proportional to ^(/ — 1)(3/ — 1), whereas 
the coupling of the 2nd order gravitational potential v''-^^ oc ^ with itself gives a kernel 

for the integrated gravitational effect proportional to "^(2/ — 1)^. This second function 
is always greater than the first one. Moreover, for quasi-flat models we expect a negligi- 
ble contribution from the coupling of lst-3rd order perturbations because / ~ 1, whereas 
for low-f2 models the integrated gravitational effect due to 2nd-2nd order perturbations is 
produced at high-z (75% of the final effect is produced in the interval [10, 10^] for = 0.1, 
see next section) where / ~ 1, and so there is no practically lst-3rd order contribution. 
At smaller z some contribution due to the lst-3rd coupling is produced but it is estimated 
to be always bounded by that due to the 2nd-2nd coupling at low-z, and this is a small 
fraction of the final contribution. 

RESULTS AND CONCLUSIONS 

We have applied the formalism of the previous section to calculate the predicted am- 
plitudes of the multipole components Ci —< af^ >. We assume an open or flat model as 
background and matter density perturbations corresponding to CDM models with either 
a primordial HZ spectrum P{k) = Ak or a LSRP one (Lyth and Stewart 1990, Ratra and 
Peebles 1994) and Qi, = 0.05, h=0.5. Primary anisotropics are normalized to the 2-year 
COBE-DMR map as given by the analysis of Cayon et al. (1995) for HZ and Gorski et al. 
(1995) for LSRP. However, secondary anisotropics are normalized to crie = 1 as this effect 
is generated by the small scale structure, ^ lOOMpc, where maybe the uie normalization 
is more appropriated (cxie is the rms density fluctuation at 16 Mpc). The results of our 
calculations are presented in the following flgures. 

We plot in Figure 1 the 2nd order power spectrum P2{k), calculated according to eq.(7), 
for different values of O and normalizing to (Tiq = 1. The amplitude of the maximum 
decreases for low Vt values, and appears at lower k. In Figure 2, we display Ci up to 
I = 1000 for different values of O, due to 2nd order anisotropics (lower curves) and due 
to primary anisotropics (upper curves). The solid, dashed and dotted lines represent the 
cases n = 1, 0.3, 0.1, respectively. For the dashed and dotted curves the results for the HZ 
spectrum are above the LSRP ones at / = 100 in all cases. With regard to the 2nd order 
anisotropics, the maximum is at / ~ 250 at the level of 10~^^ for all Q values (aie = 1). 
In the case of COBE-DMR normalization, the level of anisotropy increases by a factor 5 
for = 1 (upper solid line) and decreases by 2 and 5 orders of magnitude for 17 = 0.3 and 
0.1, respectively. 

We wish to emphasize that for Q = 1 more than 90% of the anisotropy produced by the 
second order gravitational effect is generated at a redshift z ^ 10 (Martinez-Gonzalez et 
al. 1992). The case of O < 1 is similar, for instance if O = 0.1 about 80% of the effect is 
generated at a redshift z < 30 (see Figure 3) . Since reionization could not have produced a 
substantial effect at such low redshifts, for any baryon density consistent with primordial 
nucleosynthesis, our results provide the minimal fluctuations in the CMB, independently 
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of the ionization history of the universe. Notice that for Q = Q.l there appears a knee 
in the generation of the effect with redshift which is due to the change in the expansion 
of the universe from being matter dominated to curvature dominated. The low order 
multipoles are generated during the curvature dominated phase whereas the high order 
ones are generated during the matter dominated one, as can be seen in figure (3) for 
I = 250. These two contributions at relatively high and low redshift are similar to the ones 
appearing at the linear level, the so called SW and ISW effects. 

We remark that there is no problem with the influence of non-linear scales in the second- 
order calculation because for any O the maximum of the function Ri{k) (see eq.(6)), is at 
k ^ 10. and the scales that are contributing to the multipole / have k ^ I, so for the 
calculated multipoles Ci up to Z = 1000 only scales > 6h~^ Mpc generate a 2nd order 
gravitational effect. Hence no contributions from density fluctuations with density contrast 
larger than unity enter in the perturbative calculation. 

The principal conclusion of this analysis is that a deformation {i.e. a peak) at I ^ 250 
in the 2D radiation power spectrum, C;, would be a signature of the 2nd order anisotropy 
contribution. Detection of such a feature would allow an estimate to be made of the global 
parameter Q, because the first acoustic peak emerging from the primary anisotropy is at 
1(0.) ^ 220. However, such a deformation is negligible for low-il models and it is at least 
two orders of magnitude below the corresponding primary anisotropy in the flat case. Only 
a possible experiment allowing the C^'s to be determined with an accuracy better than 1% 
could detect such secondary anisotropy. 

On the other hand, early reionization of the universe would erase the primary anisotropy 
whereas the secondary one could survive with only minor changes. Figure (4) shows 
the primary radiation spectrum for several reionization histories ranging from standard 
recombination to the case of no recombination as compared to the secondary contribution 
for n = 1 and HZ primordial spectrum. The secondary effect amounts to a contribution 
of ~ 1 — 10% in the best case. 

Finally, we remark that the second order gravitational effect is very sensitive to normal- 
ization because AT/T is proportional to the amplitude A of the primordial spectrum. 
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FIGURE CAPTIONS 



Figure 1. The second-order power spectrum for the energy-density perturbations P2{k) for il = 
1,0.3,0.1 is shown by the solid, dashed and dotted lines, respectively. The primordial 
spectrum is normalized to uie = 1. 

Figure 2. Radiation Power spectrum for the primary (upper curves) and secondary (lower curves) 
contributions. Q, = 1,0.3, 0.1 correspond to solid, dashed and dotted lines, respectively. 

The results for the HZ spectrum are above the LSRP ones at / = 100 in all cases. The 
secondary contribution is normahzed to aiQ = 1, except for the solid upper curve which 
represents the flat case with the COBE-DMR normalization. 

Figure 3. Generation of the multipole I = 250, with maximum amplitude, for the case = 
1, 0.3, 0.1 as a function of redshift (solid, dashed and dotted lines respectively). 

Figure 4. Radiation Power spectrum for several reionization models with fl — 1 and the HZ 
spectrum. Dashed lines from top to bottom correspond to the primary contribution 
in the following cases: standard recombination, r = 0.5, 1 and no-recombination. The 
solid lines represent the secondary contribution which is insensitive to the recombination 
history, upper curve is for COBE-DMR normalization and the lower one for die = 1 as 
in figure 2. 
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